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We analyze the two loop renormalization group equations in the standard model and its extensions 
for the coupling constants and the quark Yukawa couplings. The key point of our analysis is the 
observed hierarchy of the quark masses and the CKM matrix. For the one loop evolution we find 
the explicit solution for the evolution of the Yukawa couplings and show the following: 1. the CKM 
matrix depends on the energy on only one function of energy; 2. the ratios of the down quark masses 
depend on the energy through the same function as the CKM matrix; 3. the diagonalizing matrices 
of the biunitary transformation of the up quarks are energy independent. Next we give the explicit 
form of the two loop corrections to the evolution of the quark Yukawa couplings and show that the 
relative corrections are of the order A 4 for the quark Yukawa couplings and of the order A 5 for the 
CKM matrix. Finally we give the equations of the one loop evolution of the squares of the matrix 
elements of the CKM matrix and their explicit solution. 
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I. INTRODUCTION 

The problem of the Renormalization Group EVolution (RGEV) of the Quark Yukawa Couplings (QYC) and the 
Cabibbo-Kobayashi-Maskawa matrix (CKM) has been studied in many earlier papers, Refs. Despite the fact 

that many, both analytical end numerical results have been obtained on the evolution of the QYC and CKM matrix 
there is no single reference where some general properties of the evolution have been discussed. This paper is aimed 
to fill this gap and it discusses the general properties of the evolution of the CKM matrix and quark masses that 
follow from the hierarchical structure of the QYC. 

The information on the properties of the QYC come from the experimental values of the quark masses and the 
CKM matrix. Both the quark masses and the CKM matrix show the hierarchical structure with the parameter 
A« sin#c~ 0.22. This hierarchical structure must also be present in the QYC and is used to derive some conditions 
for the unknown QYC, Refs. g which reduce the number of parameters of the standard model. 

In this paper we systematically investigate the influence of the hierarchical structure on the evolution of the CKM 
matrix and quark masses by constructing the exact solution of the one loop RGEV equations compatible with the 
observed hierarchy and then considering the two loop RGEV equations with the next order corrections in A. 

The most important result that we derive is that the CKM evolution depends only on one parameter which is 
a suitable integral that depends on the model (Standard Model (SM), Minimal Supersymmetric Standard Model 
(MSSM) and Double Higgs Model (DHM)). The corrections to the one loop evolution of the CKM matrix are of 
the relative order A 5 . We next show that the evolution of the ratios of the down quark masses md/rn s and ra s /mt, 
depends on the same parameter as the CKM matrix. The evolution of the up quark masses is different. Here the 
quark masses m u , m c depend linearly on the corresponding initial values and their ratio m u /m c is constant (5) while 
the dependence for m t is non linear. 

The organization of the paper is the following. In Section II we introduce the notation and the essential experimental 
facts. Afterwards we discuss the RGEV equations up to two loops and their approximate form compatible with the 
observed hierarchy. This gives the basis for the derivation in Section III of the exact solutions of the one loop RGEV 
equations. We also present there the perturbative scheme for the study of the corrections coming from the higher 
order terms in A and the two loop contributions. In Section IV we discuss the properties of the solutions derived earlier 
and their physical importance. In Section V we derive the approximate evolution equations for the squares of the 
absolute values of the CKM matrix elements and give their explicit solution. In such a way we give the explicit form 
of the one loop renormalization group evolution of the full CKM matrix. Section VI is devoted to the conclusions. 
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II. QUARK YUKAWA COUPLINGS AND CKM MATRIX-HIERARCHY AND EVOLUTION 



The quark and lepton masses in the standard model arise through the Higgs mechanism from the Yukawa couplings 
which have the following structure in the SM 

3 

E (4 e3 ^4 + (2/u)*x<k + + h - c -)- (i) 

Here f^ e \ y u , and yd are the matrices of the Yukawa couplings of leptons and up and down quarks, respectively. <f> is 
the scalar Higgs field. The theory itself does not put any restrictions or conditions on the f^ e \ y u and yd matrices and 
from the phenomenological point of view they are constrained by the values of the physical lepton and quark masses 

Diag(m M ,m c ,ra t ) = (U u ) L y u (U u y R , Dia,g(m d ,m s ,m b ) = (U d ) L yd(U d ) R (2) 

with diagonal elements being the up and down quark masses after the spontaneous symmetry breaking and the CKM 
matrix. The diagonalizing matrices (£/ U! d) L R of the biunitary transformations transform the quark fields in Eq. (|l|) 
into the physical quark fields. As a consequence the unitary matrix (CKM matrix) 

Vckm = (Uu) L (Ud)l (3) 

appears in the charged current. 

The observed hierarchy of the quark and lepton masses and the CKM matrix is the following, Refs. (|J^] 
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with A w 0.22. 

The hypothesis of grand unification assumes that the 3 coupling constants of the standard model converge at the 
scale 10 15-16 GeV and the symmetry group of the model becomes larger (e.g. SU(5)). This may also be an origin of 
the additional symmetries or textures of the QYC at this energy scale. To relate the parameters of the QYC at the 
GU scale with the observables at low energy one conventionally uses the Renormalization Group Equations (RGE) 
for the coupling constants and the quark and lepton Yukawa couplings. The structure of the two loop RGE is the 
following Refs. 

t = j^r bl9f ~ (4^ G ' 5 ' 3 ' (6a) 



dy-u 



dt 



——8 (1) + —B {2) 



Vv,,dw (6b) 

The variable t is defined as t = ln(_E//i) and the constants bi and functions G/, [3^\ v and P^ d e v are defined for 

various models in Appendix I. The coefficients G/, /3^ d v , e u are functions of g\ and y u ,d,e,v so Eqs. (0) form a 
system of coupled non linear equations and their explicit solution is not known. The Yukawa couplings y u ,d,e,u are 
normalized in the following way 

b„]33~l, M 33 ~"A 2 , [y e ] 33 ~^~A 3 . (7) 

mt mt 

The normalization in Eq. (Q) is the origin of the hierarchy in Eqs. @ because the functions G;, (3^ d e v and (3^\ e v 
contain the squares and higher powers of y u ,d,e,v 

From Eq. (Q) it follows that positive powers oi_yd, e ,v are much smaller than the corresponding powers of y u . This 
information enables to find the solution of Eqs. (g) and to study the properties of the solutions. The other origin of 
the hierarchy in Eqs. (|^) is related to the number of loops. Here the inclusion of each additional loop is suppressed 
by the factor l/(47r) 2 ~ A 4 which is of the same order as ([yd]^) 2 and the correct approximation procedure must take 
this into account. 

Eqs. (||) will be solved and analyzed in the following steps: 
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1. Definition of the hierarchy of the equations; 

2. Exact solution of the basic equations; 

3. Perturbative corrections to the exact solution; 

4. Implications of the former results for the quark masses and the CKM matrix 

In the rest of the paper we will discuss each step in detail. 

Eqs. (|^) are the two loop RGE for coupling constants gi and Yukawa couplings y u ,d,e,u- I n these equations the terms 
of different order in A are present. As the first step we will find the approximate form of these equations neglecting 
all the terms of A and higher. These approximate equations have the following form 



bigf, (8a) 



dt (4tt) 2 

<(*) + ^y u yt + ^Tv(y u yt)]y u = j^E^y u , ( si , , 



IF = T^ K{t) + alv ^ + a 3Tr(y u2 /I)]2/ e = J^ E lVe, (8d) 
' // ' ' r Oi^(t)+a^y u yl + a^Tr(y u yl)]y v = j^-^E'(y v (8e) 



dt (4tt) 2L iV ' ZaU * U (47r) 

a . 

Eqs. (|6|) where the terms of the order A 4 are kept 



where the values of a™' d ' e ' v are given in Appendix. Eqs. (pf) will be solved exactly and as a next step we will consider 



k i 



1 rF u (yuylgi)yu = T^E?y u + -±-E%y u , (9b) 



1 (4tt) 4 ^u>^,»« (4tt) 2 (4tt) 

If = (4^^ dyd + rhf [aivdVd + a ^y^y* 

+ J^F d {y u yl 9l )y d = j^E d y d + j~^E d y d , (9c) 

+ ^(y^ 9l)Ve = J^f E tye + (i^jl^ye, (9d) 



(4^)4 www (47r)2 (47r) 

The functions F u,d ' e ' v {y u yt> 9i) are obtained from the functions (3^ d e v given in Appendix by putting yd — ye = 0. 
Eqs. (^) cannot be explicitly solved but they allow the perturbative solution by__ transforming them into integral 
equations. This allows the study of the corrections to the exact solutions of Eqs. 
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III. SOLUTION OF THE RENORMALIZATION GROUP EQUATIONS FOR y u AND y d 



The renormalization group equations (Q) will be solved in two steps. We will start by solving the one loop equa- 
tions (||). Eqs. (|8a|) are very easy to solve and their solution reads 



9l(t) 



2bigf(to)(t-t ) 



Eq. (^) is decoupled from Eqs. (|^) and ( |8cf ) and is non linear. Eqs. 
solution of Eq. (8b) is known. We therefore solve Eq. (pq) first. 



(10) 



and (|8d|) become linear once y u - the 



A. The evolution of y u from one loop RGE 



The crucial fact enabling the exact solution of Eq. (8b) is the observation that the diagonalizing matrices of the 
biunitary transformation do not depend on the energy. As mentioned before y u can be diagonalized by the biunitary 
transformation, Eq. (||). The unitary matrices (U u )l.r are the diagonalizing matrices of the hermitian matrices 
Hi = y u (t)yi(t) and H\ = y\(i)y u (t) and they fulfill the equations which follow from Eq. (pb|) 



(4tt) : 



r {aUt)+o%Hi + c4Tr(Hi)}Hi, 



1 = 1,2. 



The solution of Eq. (11) can be written in the following form 

(f - t ) k d k W u (t) 



k=l 



k\ 



dt k 



(11) 
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From Eq. (O) one can see that the derivative 



d k Hl 



dt k 



is the sum of the powers of H^(t ) with scalar coefficients. 

Thus from the hermiticity of ff*(io) it follows that H^(t) are diagonalized by the same matrices as H^(tg) so the 
diagonalizing matrices {U u )l and (U u )r of H^(t) and H^(t) are energy independent. Thus the matrix y u (t) has the 
following representation 



y u (t) = (U u ) f L A u {t)(U u ) R 



(13) 



where on the right hand side of Eq. ( [Df ) only the diagonal matrix A u (t) — Diag(m„(t), m c (t), nit{t)) does depend on t 
and {U u )l and (U u )r are constant matrices that diagonalize y u {ta)- The matrix A u (t) fulfills the following equation 
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dt " (4tt) 2 

Eq. ( |l4| ) splits into 3 equations for the up quark masses 

dm u 1 
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We first solve Eq. (15c) 
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r g (t) = exp(— - / aUr)dT). (17) 



to 



Eqs. (15a) and ( |15b|) are identical and their solution reads 



m u , c (t) = m„, c (tn)ry 2 (t)exp(^-^a 2 " / m 2 t { T )dT). (18) 



to 



The full energy evolution of y u (t) is determined from Eq. (^3|) using Eqs. ( jig ) and (jig). 

B. The evolution of yd from one loop RGE 

The one loop evolution of xu is determined from Eq. (^). The fact that the diagonalizing matrices of y u (t) do not 
depend on t simplifies Eq. ( pq) and justifies the following substitution 

Vd(t) = (^u)lW(t) (19) 

so W(t) fulfills the equation 



= (4^ {a " (<) + a " A " + "3Tr(A 2 )}!¥. 
The matrix on the right hand side of Eq. (j|^) has the following form 

{a((t)+a d 2 Al + a d 3 Tr(Al)} 



and the solution W(t) of Eq. (pp|) is the following 










/ 







: 


1 


• 


+ ^m? 





• 









\o 








Wit) = (^(t)) 1 / 2 exp(^ag /* m 2 t (r)dr) ■ Z(t) ■ W(t ) 



where 



and 



' (4tt) 2 



t 



Co 



with 



1 



(20) 



(21) 



(22) 



r' g (t) = exp(— 5 / a?(r)dr) (23) 



1 

Z(i) =(01 | (24) 




ft(t) = exp(— -^af / m 2 t (r)dr). (25) 



Putting Eqs. (|19|)-(|25|) together we obtain the following result for the one loop evolution of yd{i) 

VS) = (r;(t)) 1 / 2 (/ l (i))(^/^) (f/u) t_^ (t)(t i u)iyd(to) . (26) 
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C. Higher order corrections to the y u and ya evolution 



The two loop RGE for y u and yd, Eqs. ([)!]) and ([k]) have the following structure 

dyu,d 1 jpU.d 1 TpU.d tr,1\ 

-dT^j^r^ ^+(4^^ y - d - (27) 

In Eqs. ( p7| ) the terms containing the functions are suppressed by the factor l/(4ir) 2 ~ A 4 in comparison to the 
terms with E'"' d . Eqs. @> can be transformed into the integral equations 

W»,«i(t) = VuAto) + ^ jf ^'V^)^ + ^ jf E^ d y u4 {r)dT (28) 

and Eqs. ( pq ) can be transformed into the following recurrence relation 

yJ2(t) = lfa.d(*d) + ^ £ ^W^Wdr + ^ £ ^^(rjdr. (29) 

The functions y^\{t) converge for n — > oo to the solution of Eq. (|27|) for an arbitrary initial function y^\{t). If we 
choose y^\{t) to be the solution of the one loop RGE, Eqs. ( |l3|) and (|26| ) then we obtain 

»£!(*) = fa.-(*o) + ^ f K'y^r + j±y 4 f E? d y^dT. (30) 



to 



Now 



lM(«o) + ^ / ^'"l/S^ = ^(*) (31) 



because and fulfill Eqs. (8b) and (|8c[), respectively. It thus follows that 



»S(*)=A*) + 7Al /V*/i> (32) 



(^) 4 y* 

which means that the lowest order corrections to the solution of the one loop RGE for y^\{t) are of the relative order 
A 4 when compared to the y^\{t) and Eq. (|32|) gives the explicit form of this correction. 



IV. PROPERTIES OF THE RENORMALIZATION GROUP EVOLUTION OF THE QUARK MASSES 

AND THE CKM MATRIX 

In this section we will study the physical implications of the results obtained in the previous section. We will 
present them in the series of the theorems 

Theorem 1 The one loop RGEV of the CKM matrix depends on only one function of the energy h(t), given in 

Proof: The CKM matrix, given in Eq. ||) is constructed from the matrices (U U ) L and (Ud) L - In section III. A 
we have shown that (U U ) L does not depend on the energy. The dependence on t can only be contained in (Ud) L 
which is determined from the matrix yd(t) of the down QYC. The matrices {Ud) L and (Ud) R do not depend on the 



normalization of yd(t). From Eq. (26) one can see that the only dependence of yd(t) on t apart from the normalization 
is contained in the matrix Z(t), wl 



rhich is the diagonal matrix and as can be seen from Eq. ( |24| ) it depends only on the 
function hit). This completes the proof of Theorem 1. 

Theorem 2 The ratios of the down quark masses md/rn s and m s /mb are the functions of only h(t), given in Eq. (fe 



G 



Proof: The down quark masses are obtained after the diagonalization of the down QYC y d (t) and the ratios of m,d/m a 
and m s /mb do not depend on the normalization of y d (t). It thus follows from the proof of Theorem 1 that these ratios 
are the functions of h(t). 

Theorem 3 The ratio of the up quark masses m u /m c is energy independent. 
Proof: This theorem follows directly from Eq. (|l^) . 

Theorem 4 The next order correction^of the RGEV of the CKM matrix are of the order A 5 . 

Proof: If we consider RGE in the next order then the terms of the order A 4 are preserved in the equations and the 
evolution is governed by Eqs. (||) and the explicit solution in the next leading order is given by Eq. ([32]). From Eq. ( |32] ) 
one obtains by the direct calculation the following result for the commutators 

[ftjMVA 5 , (33a) 

[y ( Sv%ySvZ]-^ (33b) 

From Eqs. (|3^) using the time independent perturbation theory one obtains that the diagonalizing matrices of the 
biunitary transformations in the next order are corrected by the terms of the order A 5 and this completes the proof. 

V. EVOLUTION OF THE CKM MATRIX 

In this section we will find the equations for the one loop evolution of the squares of the absolute values of the 
off-diagonal elements of the CKM matrix. 

The CKM matrix is defined in Eq. (H). Using the fact that {U u )l is energy independent we obtain from Eq. 12( 



(U u ) L y d (t)y d (t)Hu u )[ = r;(t)(/ l (t))( 2 "3M)z(t)(C/ u ) i?/d (t )y rf (t )t(C/ u )[z(<) (34) 
which can be written in the following way 

VcK M {t)M 2 d {t)V^ KM {t) = r;W(/i(i)) (2a 3/^)z(i)(C/„) L y £i (i ) 2/(i (t )t(C/ u )[z(t) (35) 

where Mj (t) is the diagonal matrix of the squares of the physical down quarks Yukawa couplings which become the 
squares of the down quark masses after the spontaneous symmetry breaking. Now differentiating Eq. ( |35| ) with respect 
to t we obtain the following result 

vX (t) dVcKM 

= {M 2 d )- x V^ CKM {t)^^Ml - (Mj)- 1 ^ ifM (t)|((C/ u ) i y rf (i)y d (t)t(C/ u )t / )y CKM ( t ) + {Mjy^. (36) 



The second term on the right hand side of Eq. (|36[) is equal 



(37) 



where the vector R = (Vtd, Vt s , Vtb) and h(t) is given in Eq. (|2 

The off diagonal matrix elements of VQ KM (t) dVc d f M can now be evaluated from Eq. (|36"1). Using these matrix 
elements we obtain the following evolution equations 

d \ V "b\ 2 _ h' mg + mg fo' mg+rog _ 2h' , |2iy 12 

T7 — _ T 2 2y V ub V udV td V tb ) - —— 2\ V ub V usV ts V t b) + C.C T-\V u b\ \V t b\ , (38aj 

at h m d — ti mj. — h 



1- we include in RGEV equations the terms of the order A 4 that come from the one loop RGE that contain the terms y d y d and 
the two loop contributions of the order ~ 1 multiplied by l/(47r) 4 . 
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d\V r h\ h mi + mr, „ h mi + ml , „ ^ 2/i , , , , , 

= -Tl^t {v '^ v ^ - ¥sM (n,li - V " V,t) + " * -X^l |v »l • (38b) 

d \ V cd? _ ti ml + m\ h' ml + m\ 2h' 2 2 

-^-~T^f^ ( ^ Ks ^ d) ~¥;^^ |Ftd| ' (38d) 

Eqs. (|38| ) can be solved explicitly and the solution reads 



W = i T/0 i 2 ; L f 7VTT . ( 39b ) 



r^°, | 2 

1 liD 1 


l^| 2 (^ 2 -l)- 


f 1' 






|vSl a C^ 2 -i)- 


f 1' 






|V&| 2 (l-fc a )H 


-/i 2 






|V&| 2 (l-fc 2 H 


-/i 2 



l^l 2 = iTTcT^TirTT^ (39c) 



|T/cd|2 " IOTi-V> + /> 2 ' (39d) 

Here |V^-| 2 are the initial values of the squares of the absolute values of the corresponding CKM matrix elements and 
h is given in Eq. (^5|). From the unitarity of the CKM matrix and using Eqs. (j3^) one can calculate the absolute 
values of all the remaining matrix elements of the CKM matrix and thus we determine the renormalization group 
evolution of the full CKM matrix. 

VI. CONCLUSIONS 

In this paper we have analyzed the renormalization group evolution of the quark Yukawa couplings and the CKM 
matrix based on the observed hierarchy of the quark masses and the CKM matrix. The inclusion of the hierarchy 
greatly simplifies the analysis and leads to simple, explicit results for the evolution of the Yukawa couplings (Eqs. ( |l3| ) 
and Eqs. (|26|)) and the CKM matrix (Eqs. (p9|)). The other remarkable result is that the diagonalizing matrices of 
the up quark Yukawa couplings are energy independent in the leading order. This means that the transformation 
{^u)l,r — > (U u )l,r(iPu)l,r, {i>d)L,R — > (U u )l,r(iP(i)l,r wi U diagonalize the matrix of the up quark Yukawa couplings 
and it will stay diagonal upon the renormalization group evolution and the CKM matrix will be determined only from 
the down quarks Yukawa couplings. This fact may simplify the model building based on the symmetries of the quark 
Yukawa couplings. 
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APPENDIX: CONSTANTS IN THE RENORMALIZATION GROUP EQUATIONS 



The two loop renormalization group equations for various models have the following structure 
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where 6/ are constants dependent on the model and G; , and are the following functions of the coupling 
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The values of the coefficients bi and a l k for the one loop renormalization group equations are given in Tables [j], || 
and III (see e.g. pi- 
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TABLE II. The coefficients a{ for various models. The 
model dependent constants a, b and c are given in Table til. 
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TABLE III. The constants a, b and c for various models. 



Model 
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SM 


1 
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-1 
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The coefficients for the two loop renormalization group equations that appear in the functions Gi and j3\ can be 
found in Ref. §. 
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